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entangled number state may be performed. They lead to a violation of a Clauser-Horne Bell inequality for 
any finite photon number N. We also present an analytical expression for the two-mode Wigner function of a 
maximally path-entangled number state and investigate a Clauser-Horne-Shimony-Holt Bell inequality for such 
a state. We test other Bell-type inequalities. Some are violated by a constant amount for any N. 
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I. INTRODUCTION 



Maximally path-entangled number states of the form 
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V2 



(\N)a\0) b + J*\0)a\N) b )> 
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(often referred as NOON states) have important applications to 
quantum imaging (J], metrology 10, 13], and sensing [4]. Char- 
acterizing their quantum mechanical properties is therefore a 
valuable task for improving upon the above applications. En- 
tanglement is the most profound property of quantum me- 
chanical systems. NOON states are non-separable states and 
hence are entangled. But do they also show nonlocal behav- 
ior when we perform a correlation experiment on their modes? 
The amount of nonlocality demonstrated by a Bell-type exper- 
iment provides an operational definition of entanglement (for 
a review of Bell inequalities and experiments see, e.g., 051]. It 
distinguishes between the class of states that are entangled but 
admit a local hidden variable model and those which do not 
and so may be called EPR correlated |6|] . 

Several publications [7] address the question of whether the 
NOON states are EPR correlated for the case N=l. Gisin and 
Peres have shown that it is possible to find pairs of observ- 
ables, whose correlations violate a Bell's inequality for any 
nonfactorable pure state of two quantum systems ^ • This re- 
sult was later extended to states of more than two systems by 
Popescu and Rohrlich |0] . Recent experiments [03 EH] have 
reported strong evidence that NOON states violate a Bell's in- 
equality for N = 1 , leaving open the question as to what ex- 
periments might show EPR correlations for N > 1 . 

We propose a specific experiment that shows that NOON 
states are EPR correlated for any finite N. We investi- 
gate two measurement schemes using the unbalanced homo- 
dyne detection scheme described in Ill2ll and compare the re- 
sults. The correlation functions we calculate can be related 
to two well-known phase space distributions: the two-mode 
Q function and the two-mode Wigner function. Banaszek and 
Wodkiewicz first pointed out the operational definition of the 



Q and Wigner function 11211 . We modify this approach and 
calculate the distribution functions for the NOON states en- 
tirely from these phase space distributions, and thereby con- 
struct a Clauser-Horne and a Clauser-Horne- Shimony-Holt 
Bell inequality. In section [TV] we also test other Bell-type in- 
equalities not commonly used so far in quantum optical ex- 
periments. 

The above Bell-tests may be performed in an unbal- 
anced homodyne detection scheme as given, for example, in 
Ref. lfl2h and shown in Fig.Q] For simplicity we choose (p = 71 
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FIG. 1: Unbalanced homodyne detection scheme for a Bell experi- 



ment with NOON states. Here I 1 ?) 
and b label the modes. 
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for the states in Eq. (Q]). 

It is now understood that the introduction of a reference 
frame is required in any Bell test fl3ll and one should con- 
sider the field modes as entangled rather than the photons 
In the number basis, a shared local oscillator acts 
as the required reference frame. The beam splitters in this 
approach are assumed to operate in the limit where the trans- 
mittivity T — > 1. We further assume that a strong coherent 
state |y), where |y| — > °o, is incident on one of the two input 
ports. The beam splitter then acts as the dis plac ement operator 
D(yy/1 - T) on the second input p ort iflej \vK [l8ll . We intro- 
duce complex parameters a = J a V^ — T and (3 = y^Vl — T. 
The phase- space parameterization with respect to these pa- 
rameters is then analogous to a correlation experiment with 
polarized light and different relative polarizer settings where 
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the nonlocality of polarization entangled states such as = 
(\H) a \V) b - \V) a \H) b )/y/2is well established. 



II. BELL EXPERIMENT WITH ON-OFF DETECTION 
SCHEME 

In the first experimental setup we consider a simple non- 
number resolving photon-detection scheme. In the case of 
the homodyne detection scheme under consideration, the lo- 
cal positive operator valued measure (POVM) is defined by 
Q(a)+P(a) = i with, 



G(a) = D(a)|0)(0|Dt(a), 

oo 

P(a) = D(a)^\n)(n\D\a). 



(2) 
(3) 



n=\ 



We assume lossless detectors for our investigation. The ex- 
pectation value of 2(a) tells us the probability that no pho- 
tons are present, depending on the phase and amplitude of 
the local oscillator. The expectation value of P(a) gives the 
probability of counting one or more photons, while not dis- 
tinguishing between one or more photons. So we simply as- 
sign a 1 to a detector click and a otherwise, giving us a 
binary result. We label the two modes of the NOON state by 
a and b. The corresponding measurement operators for a cor- 
related measurement of the displaced vacuum can be written 
as G fl ((x) Gfc(P)- The expectation value for the state is 
given by 



G^(a,P) = <¥|G fl (a)<8>G*(P)pF> = |<a,P|¥)| 2 



(4) 



The above expression is the two-mode Q function of the NOON 
state up to a factor l/ft 2 , and the result is given by 



= ^- (|a|2+lp|2) |a /V -P /v | 2 . (5) 



To obtain the probabilities for the individual measurements 
we calculate 

Q a (a) = (V\Q a (a)®i b \V) = ^e-M 2 (M^ + i),(6) 

g*(P) = (»P|i fl ®a(P)W = ^- |p|2 (^ + i).(7) 

Using the completeness relation P(a) = 1 — G(°0> we obtain 
the probabilities for the correlated and single detector counts 
— P a (a) = 1 - G fl (a), ft (P) = 1 - G*(P), and ft*(a, p) = 1 - 
Ga(oc) - Gfc(P) + Gafe (ot, P) — in terms of the Q functions. We 
build from these the Clauser-Horne combination (CH) lfl9ll . 
which for a local hidden variable model admits the inequality, 



-1 < ft^(a,P)-ft^(a,p / )+ft^(a / ,P)+ft^(a / ,p / ) 
-ft(a')-ft(p)<0. 



(8) 



If this inequality is violated it follows that NOON states contain 
nonlocal correlations, i.e., EPR correlations. In order to attain 
such a violation, we minimize the function CH = ft&(oc, P) — 



Paft(a,P') + ft*(a',P) +ft fe (a / ,p / ) -ft (a') -ft(P) for a 
given N over the parameter space spanned by a, a 7 , p, and 
P 7 . The violation of the Clauser-Horne combination for the 
NOON states with N = 1 , . . . , 4 is shown in Fig. [2] The results 
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FIG. 2: Violation of the Clauser-Horne Bell in- 
equality as a function of N. 

show a decrease in the amount of violation with N. The max- 
imal violation is obtained for N = 1 . For N > 3 the violation 
is so reduced that it would be increasingly hard to observe ex- 
perimentally. If we increase the precision of our numerical 
method, we observe that for large N the minimum of the CH 
combination, in fact, never hits the classical bound of — 1 ex- 
actly, i.e., there is a violation of the inequality for any finite 
N, which can be shown as follows. Let N be finite and odd. 
We choose a! = P = and a = — p 7 , then the CH combination 
reduces to CH = \/N\ |a| 2A ^H a l 2 (l - 2e~^ 2 ) - 1. For any 
< |a| 2 < ln2, we obtain CH < -1. For even N the same 
proof holds except that we need to choose a = p 7 instead. 

The Bell measurement presented leads to a decrease of the 
amount of violation with N. This decrease with N is due to 
the specific way the reference frame is introduced in terms 
of the local displacement operators D(a) and D($) for the 
correlation measurement. The scheme is based on measuring 
the overlap of coherent states with the modes of the NOON 
state. The elements contained in Eq. © are of the form 
(Af,0|oc,P) and (0,Af|oc, p). In order to maximize those prod- 
ucts we would need a to take, at the same time, the values 
|a| 2 = N and |a| 2 = 0. Since the 'distance' of N to the vac- 
uum becomes larger with N, the correlated overlap is reduced. 
This may explain the decrease in the amount of violation ob- 
served. 

We can also display some correlations by plotting the 
marginals of the Q function in Eq. Q. We therefore decom- 
pose the dimensionless complex local oscillator amplitudes in 
the set of real variables x, y, w, v, i.e. , a = x + i y, P = u + i v, and 
obtain Q m (y, v) = f™^ f™^ Q a b (x, y,u,v) dxdu. These proba- 
bility densities are displayed in Fig. [3 01 and[5l for N = 1 , 2, 3. 
We see that the distributions for N = 2, 3 have a higher sym- 
metry than for N = 1 . 

The linear correlation coefficient r = cov(y,v)/(AyAv), 
where cov(y, v) = f^y -y)(v- v)Q m (y, v) dy dv, van- 
ishes for all N > 1 , although we see from the pictures that the 
two phase space variables are statistically dependent. This is 
an indication of nonlinear correlations between the two phase 
space variables. Note that the measurement described by the 
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FIG. 3: The marginal Q function Q m (y,v) forN= 1. 
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FIG. 4: The marginal 2 function Q m (y, v) for N = 2. 



operators in Eq. and Eq. © requires only non-number re- 
solving photon counters and may therefore be performed with 
current detector technology. In the next section we consider 
a correlated parity measurement on the modes and investigate 
the amount of violation in this scheme. 



III. BELL TEST WITH PARITY MEASUREMENT 

An operational definition of the two-mode Wigner function 
for the NOON state is given in terms of a correlated parity mea- 
surement [12 1. The measurements can be described by the 
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FIG. 5: The marginal Q function Q m {y, v) for N = 3. 



following POVM operators: 

oo 

ft+(a) = fi(a)J^\2k)(2k\ff(a), (9) 

k=0 

ft" (a) = D(a) £ |2^+l)(2^+l|D t (a). (10) 

k=0 

The corresponding operator for the correlated measurement of 
the parity on mode a and b may be defined as: 

n(a, p) = (ni +) (a) - fijf >(a)) ® (ni +) (P) - ft£->(P)) . 

The outcome of the measurements are either + 1 or — 1 . It may 
be noted that this operator can be rewritten as 

n(a,(3) =D a (a)D b (m-lf a+ " b Dl(a)Dl^), (11) 

and is equivalent to the operator for the Wigner function in 
HH (up to a factor 4/% 2 ). We note that the operator in 
Eq. (ITTI) is essentially a product of operators for mode a and 
b: 

ft(a,|3) =A fl (a)(-l)^J(a)^(P)(-l)^(p). (12) 

Using this property the expectation value of Eq. (fT2l) for the 
NOON state can be expressed as a function of two Laguerre 
polynomials and an interference term, 

n(a,p) = (v|n(a,p)|v) = 

|^- 2 l-l"- 2 IPI 2 [(-l)^(^(4|oc| 2 ) -r-^(4|P| 2 )) 

-±-(a* N V N + a N p* N )}, (13) 

where Ln(x) is the Laguerre polynomial i22ll . The two-mode 
Wigner function is obtained from W(a, P) = n(a, P)4/7i 2 . By 
building the CHSH Q inequality defined by 

- 2 < n(cc, P) + n(a 7 , P) + n(a, PO - n(a 7 , PO < 2 , (14) 

we determine how this Bell inequality is violated as a function 
of N. A minimization procedure in the parameter space a, a 7 , 
P, and P' as a function of N is carried out with a numerical 
routine to investigate the amount of violation. We see that 
the correlated parity measurement leads to a violation of the 
CHSH Bell inequality for N = 1, and that states with larger N 
do not violate the inequality. 

The Wigner function may also be used to understand this 
behavior. We therefore calculate the marginals of the Wigner 
function by integrating over two of the variables, where we 
use the same decomposition of the dimensionless complex lo- 
cal oscillator amplitudes a = x+iy and a = u+iv, and obtain 
Wm(.y, v) = /Too f-ooW(x,y, u, v)dxdu. The function W m (y, v) 
is positive definite and can be interpreted as the probability 
density for the remaining variables. From the density plots in 
Fig- IS and [8] we see that the probability densities become 
more symmetric the larger N becomes, similar to the previous 
case for the marginals of the Q function, but the interference 
structures are much more pronounced than for the Q function. 
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FIG. 6: The marginal Wigner function W m (y, v) for N = 1. 
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FIG. 7: The marginal Wigner function W m (y, v) for N = 2. 



Here we also obtain a vanishing correlation coefficient r for all 
N > I, from which we can infer that a nonlinear correlation 
measure is necessary to describe these correlations. 

We conclude from the results of the first section that a set 
of parameters can always be found which violate the CH in- 
equality in Eq. ([5]). Therefore NOON states show EPR correla- 
tions for any finite N. The presented setup is not yet optimal 
but might be promising for demonstrating EPR correlations of 
NOON states with low photon numbers N experimentally. Al- 
though, the requirements for the overall detection efficiency 
for a loophole-free test of the CHSH Bell inequality would be 
very large, i.e., 96% for N = 1 ifTTll. 

In the following section we are going to show that the test 
of other Bell-type inequalities leads to a different result. 
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FIG. 8: The marginal Wigner function W m (y,v) for N = 3. 



IV. MORE BELL-TYPE INEQUALITIES 

So far we have used the CH and the CHSH Bell inequalities 
defined in Eq. ^ and Eq. (IT4l> . Other Bell inequalities might 
be more suitable for a Bell test for a nonlocal experiment with 
NOON states. The CH Bell inequality is a specific inequality 
for four correlated events, where at most two are intersected 
at the same time. Pitowsky [24] derived all the Bell-type in- 
equalities for three and four correlated events: 

0<Pi- Pij ~ Pik + Pjk, (15) 
Pi + Pj + Pk-Pij-Pik-Pjk < 1, (16) 

-l<Pik~P jt + Pii + Pjk~Pi-Pk<0, (17) 

for any different ij,k,£. Eq. (TTtT) is the CH inequality. 
Eqs. (115116b are inequalities in the so called Bell- Wigner poly- 
tope of three correlated events, whereas Eq. (fTTl) belongs to 
the Clauser-Horne polytope ll24ll . Later on Janssens et al. ll25ll 
explicitly constructed inequalities for six correlated events 
where, as before, two are intersected at the same time. We 
consider the following four: 

Pi + Pi + Pk + Pt~ Pij ~ Pik~ Pit- Pjk- Pj£~ Pkt < 1, (18) 

2pi + 2pj + 2p k + 2pi - p^ - p ik - pa -p jk -pj£- pu < 3, 

(19) 

< Pi - Pij - Pik ~ Pit + Pjk + P jt + Pkt , (20) 

Pi+Pj+Pk-2pt-Pij-Pik+Pit-Pjk+Pjt+Pkt < I- (21) 

We investigate the amount of violation for the inequalities in 
Eqs. (1181121b for the simple on-off detection scheme of section 
[J with the detection probabilities given by Eqs. (1516171) . The 
probabilities in Eq. dTSt are then replaced by 

Ji = 2(a) + g(P) + G(7) + 2(8) - G(«, P) - G(a, y) 

- e(a,8) - G(P,y) - G(P,5) - G(y,5), (22) 

so that the inequality is given by J\ < 1 . We make the follow- 
ing assignment i — > a, j r — > P, k — > y, and t — > 8. The single- 
count probabilities G(°0 can either be measured by Alice or 
by Bob. The joint probabilities are always measured between 
Alice and Bob. 

A maximization procedure carried out in the parameter 
space a, P,y,8 leads to a constant violation of the inequality 
as shown in Fig. 13 This new result will be interpreted in more 
detail at the end of this section together with the results from 
the remaining inequalities. 

The probabilities in Eq. (IT9l) can be rewritten in terms of 
the local oscillator amplitudes as well 

h = 2Q(a) + 2G(P) + 2Q(y) + 2Q(8) - G(a, P) 

- G(a,y) - G(a, 8) - g(P,y) - G(P, 6) - G(y, 8) ,(23) 
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FIG. 9: Violation of the inequality Eq. (l22t as a 
function of N. 
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FIG. 1 1 : Violation of the inequality Eq. (|25l > as 
a function of N. 



where the inequality is then given by J 2 < 3. A maximiza- 
tion procedure for the parameters in Eq. (l23t shows a constant 
violation of 4 for any N. 

Finally the probabilities in Eq. (l20l) and Eq. (l2Tb appear to 
be, in terms of the complex parameters a, (3,y, 8, 



on her/his particle as given by 



h = G(a) ■ 



■G(«,P)- 
-G(P,Y)- 



■G(a,Y)- 
-Q(P,8)- 



-Q(a,8) 
■G(Y,8), 



(24) 



with the inequality < J3 . And 



/ 4 = 2(a) + G(P) + G(Y) - 22(8) - Q(a, (3) 

- G(a, Y) + G(a, 8) - g(P, Y) + G(P, 8) + G(Y, 8) ,(25) 

with the inequality J4 < 1 . Unlike the two previous cases we 
do not obtain a constant violation for Eq. (l24l) . Instead we 
attain a decreasing violation with the photon number N as dis- 
played in Fig.[lOj So not all inequalities in the polytope of six 
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FIG. 10: Violation of the inequality Eq. (1241) as 
a function of N. 

correlated events can be violated by a constant amount. How- 
ever the last inequality Eq. (l25t is violated constantly again 
with a value of 1.5 as displayed in Fig.fTTI 

The Bell-type inequalities with six correlated events all 
show a stronger violation than the CH and CHSH inequali- 
ties. We attain, except for one case, a constant violation for 
any N. 

We expect that Bell inequalities exist which show a constant 
violation because of the following argument. Let's assume 
Alice and Bob can perform locally, a unitary transformation 



Uj = \l)j{N\ + \N)j(l\+ £ \n)j(n\ 

n=0 



(26) 



where j = a,b. The combined application of their local 
unitary transformations transforms the one-photon entangled- 
state into an N-photon entangled NOON state 



U a U b 



1 1,0> — |0, 1> _ \N,0) - \0,N) 



(27) 



V2 V2 
(see also Fig. [12]). The fact that this local unitary operation 



, D a {a) 



ihiff) 



|7«) 
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FIG. 12: Alice and Bob apply a local unitary operation on her/his 
mode a and b. 

exists tells us that there ought to be a nonlocal measurement 
which acknowledges this fact. Therefore the same amount 
of nonlocality should be obtained for the Af-photon state as 
for the one-photon entangled state. The fact that some of the 
Bell-tests do not show this result means that these Bell-tests 
are not optimal. However the Bell-tests of the inequalities in 
Eqs. (122123125b seem to be optimal for the NOON state since 
their outcome shows a constant violation for any N. We point 
out that these Bell-type inequalities have two more joint prob- 
abilities than the CH and CHSH Bell inequalities. The class of 
inequalities with six joint probabilities seem to be more sen- 
sitive to the nonlocality in NOON states. From our results we 
also infer, that for some applications, types of Bell inequal- 
ities other than the Clauser-Horne and the Clauser-Horne- 
Shimony-Holt should be considered. It is, experimentally, not 
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more difficult to test these Bell inequalities; since one only 
needs to measure the correlation functions for a few more pa- 
rameter settings. The experimental setup does not need to be 
changed. 

V. CONCLUSION 

We presented several Bell-tests for NOON states. In sec- 
tion [III a simple on-off detection scheme together with the CH 
Bell inequality shows a violation for any N although the vio- 
lation decreases as N increases. In section Hill we consider a 
correlated parity measurement together with the CHSH Bell 
inequality. A violation is found only for N = 1 . In section 
HV1 we consider the simple on-off detection scheme but test 
Bell-type inequalities with six joint probabilities. We then at- 
tain a violation that stays constant for any N and we show 
by a simple argument with local unitary operations that this 
is to be expected for an optimal Bell-test with NOON states. 
If we use the violation of a Bell-type inequality as a measure 
of nonlocality then NOON states contain the same amount of 
nonlocality for any N. Despite this fact, using NOON states 



with large N is advantageous for applications like quantum 
imaging, metrology, and sensing, although the improvement 
in the performance of these applications does not seem to be 
necessarily related to the nonlocal properties of NOON states. 

Finally we point out that it might be advantageous in many 
experiments to also test the Bell-type inequalities in section 
HVl in addition to the CH or CHSH Bell inequalities. One 
gains more insight into the nonlocal properties of the states 
under investigation, as shown by our example. 
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